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O ' Abstract 



We study properties of the oscillation effects in the matter of the Earth on antineu- 
trino fluxes from supernovae. We show that these effects can provide explanation 
of the difference in the energy spectra of the events detected by Kamiokande-2 and 

,^ . 1MB detectors from SN1987A as well as the absence of high-energy events with 

E ^ AO MeV. This explanation requires the neutrino oscillation parameters Am? 
and sin^ 20 to be in the region of the LMA solution of the solar neutrino problem 
and the normal mass hierarchy if |C/e3p ^ 10"'^. The hierarchy can be inverted if 

5^ \ |f^e3p ^ 10~^. The solution of the solar neutrino problem based on fe-conversion 

to a pure sterile state is disfavoured by SN1987A data. 



1 Introduction 

The detection of the neutrino burst from the supernova SN1987A by the Kaniiokande-2 
|l|] and 1MB [Q] detectors has confirmed the general picture of gravitational collapse, hot 
neutron star formation and neutrino emission (see p, ^ for a review). 

At the same time certain features of the detected neutrino signals remain unexplained. 
In this paper we will consider the difference in the energy spectra of the events detected 
by Kamiokande-2 (K2) and 1MB. It is difficult to explain such a difference in terms of the 
characteristics of the detectors (energy thresholds, efficiencies of detection, sizes, etc.). 
The i>e spectrum inferred from the K2 events is substantially softer than that from the 
1MB events. This can be quantified by comparing the parameters of the original Ve 
spectra, such as the effective temperature Tg and the luminosity Lg, that were estimated 
from the results of the two experiments. It was found that the regions in the Tg — Lg plane 
determined from the K2 and 1MB data have only marginal overlap and the probability 
that the two sets of data correspond to the same spectrum is less than few per cents [^, |^ . 

More specifically, we will discuss the following features of the SN1987A signals: 

(i) Concentration of the 1MB events in the energy interval i? ~ 35 -f- 40 MeV. 

(ii) Absence of events at 1MB above E' ~ 40 MeV (which looks like a sharp cut of the 
spectrum) . 

(iii) Absence of events with i? ^ 35 MeV at K2. 

The effects of mixing on supernova neutrinos have been extensively studied [0, ||. In 
particular, soon after the observation of SN1987A it was marked that the differences in 
the K2 and 1MB spectra could be related to oscillations of v^ in the matter of the Earth 



and to the different positions of the detectors at the time of detection p|, |T^. It was 
realized that, for this oscillations mechanism to work, one needs lS.w? ~ 10~^ eV^ (i.e. in 
the region of the Earth regeneration effect) and large (close to maximal) mixing of Ve- 

At that time the idea did not attract much attention, since a large lepton mixing 

was considered unnatural and conversion with small mixing angle was indicated as the 

most plausible solution of the solar neutrino problem. Moreover, later it was argued 

TT|, O that a large mixing of electron neutrinos is disfavoured by the SN1987A data: the 



^e <-*■ Vy^jvr conversion inside the star leads to the appearance of a high energy tail in the 
spectrum which contradicts observations. The bound is absent, however, in the range of 
Earth matter effects, lS.rr? ~ 10"^ eV^. 



Later, the likelihood analysis of the combined K2 and 1MB spectra in ||T3[ has con- 
firmed the results of |T^. In general, the fiavour conversion of v^ leads to worser agreement 
with the supernova theory. Some mild improvement of the likelihood appears for oscilla- 
tion parameters in the Large Mixing Angle (LMA) region when the Earth matter effect 
is taken into account. This improvement is obtained, however, at the price of aggravat- 
ing the confiict with the SN theory: it was found that the best fit point corresponds 



to too low temperature of the original z/g spectrum, Tg ^ 1.9 MeV, and too high total 
binding energy, L ~ 10^^ ergs. Moreover, the best fit values Ath? ~ (3 -^ 5) ■ 10~^ eV^, 



found in |13| are now excluded by the Superkamiokande data on the day-night asymmetry. 



Since the first proposal [^j the situation has drastically changed. With very high 
confidence level the atmospheric neutrino data are explained by v^ ^^ v^ oscillations with 



maximal or nearly maximal mixing |]14|. The LMA solution gives the best global fit of 
all the available solar neutrino data |15|. It seems that large mixing is a general property 
of leptons. Clearly, if the LMA is the true solution of the solar neutrino problem, then a 
significant part of the z/g events detected from SN1987A were produced by the converted 
muon and tau antineutrinos. This means that in 1987 we observed the first appearance 
signal of neutrino conversion! 

In view of this we have reconsidered the explanation of the features of the neutrino 
signals from SN1987A in terms of neutrino conversion in the star and in the Earth. 

In section ^ we study the general properties of the conversion of antineutrinos in the 
matter of the supernova and of the Earth. We apply the results to the explanation of 
the difference of signals detected by Kamiokande-2 and 1MB in section 0. In section ^ we 
summarize implications of the Earth matter effect on antineutrinos for the neutrino mass 
spectrum. 

2 Neutrino conversion in the star and in the Earth 

We assume that the electron neutrino mixes with some combination of the muon and 
tau neutrinos, z/^, (z^e = cos6'z/i + sin6'z/2, v'^ = cos6'z/2 — sin6'z/i) and the oscillation 
parameters of the Ug — u' system are in the region of the large mixing angle solution of 
the solar neutrino problem: 



sin^ 26 ~ 0.6 ^ 1.0 , Am^ ~ (1.5 H- 10) • 10~^ eV' 
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(the case of the LOW solution will be considered later). We will also assume that other 
Am^ and mixings in the neutrino spectrum are irrelevant for the antineutrino channel. 
With the parameters ([^) the propagation of antineutrinos inside the star is completely 



adiabatic, so that the following transitions occur |jT6|: 



Z/e ^ Z/i , Z/' -^ Z/2 . (2) 



That is, the originally produced electron antineutrino will reach the surface of the star, 
and consequently the surface of the Earth, as the pure mass eigenstate z/i, and u'^ as the 
pure mass eigenstate z/2. As a result, the flux of the electron antineutrinos at the surface 
of the Earth equals: 

Fg = cos" OF^ + sin^ ^F° = F° - sin^ e{F^ - F° ) , (3) 

where Fg and F£ are the original fluxes of z/g and u' produced in the center of the star. 



In the matter of the Earth the mass eigenstates ui and z/2 oscillate. Taking into account 



these oscillations we find the fiux of i/^ at the detector [0 : 



Ff = F.+ (Pie--cos2^)(F°-F°), (4) 

where Pig is the ui -^ z/e transition probability in the matter of the Earth; in the derivation 
of eq. (^ we used also the relation P2e = 1 — Pie- Since for both K2 and 1MB the neutrino 
trajectories were in the mantle, where the density change is rather small, we use for Pig 
the formula for oscillations in uniform medium [1121 : 



Pig = cos^ 9 + sin 2em sin {29 - 29m) sin^ ( 7- I • (5) 



'"in 



Here d is the length of the neutrino trajectory inside the Earth for a given detector, 9m 
and Im are the mixing angle and the oscillation length for the antineutrinos in the matter 
of the Earth: 



.29/1 sin 2g _ sin2^ 

'''' "^ {cos29 + {l/loW + sm^29' " " sm29 ' ^^ 

In eq. (H) l^ = AnE/Am? and Iq = 2ij: / {\/2G pne) are the vacuum oscillation length and 
the refraction length; Gp is the Fermi coupling constant and rif. is the electron number 
density. In the expression (||) 9m and Im should be taken for the average electron number 
density rig along the neutrino trajectory for a given detector. 

The qualitative features of the conversion effects can be immediately seen from eqs. 
(0)"®5 taking into account that the original v' flux has harder spectrum than the one of 
Ve and the total luminosities in the different flavours are comparable. This means that a 
critical energy, Ec-, exists such that: 

F^, > F9 for E> Ec 

P^°, < F^ for E <Ec. (7) 

The energy E^ depends, in particular, on the effective temperatures and luminosities of 
the original u' and z/g fluxes: E^ = Ec{Tg, T^ii, Lg, L^')- 

According to eq. (0) the conversion in the star leads to a composite spectrum of 
electron antineutrinos with sin^ 9 admixture of the hard component which generates an 
high energy tail in the sample of events. In comparison with the original u^ spectrum. Eg, 
the flux Pg is suppressed at E < Ec and enhanced at P > Ec. The effect is proportional 
to sin^ 9 and the conversion probability itself does not depend on the neutrino energy. 

The oscillations inside the Earth give an opposite effect (see eqs. (§), (|)): they 
enhance the flux at E < Ec and suppress it at high energies. Notice that, since Pig > 
cos^ 9, for LMA parameters and in constant density approximation the matter of the Earth 
always regenerates the Uf,'- the oscillatory factor in eq. (^ is positive since 29m = (0 H- tt) 
and 9 > 9m- Shortly, an increase of the medium density will tend to return the neutrino 



state into the initial state, that is, z/g at very high densities. Entering the Earth is 
equivalent to such an increase. 

Combining eqs. (Q) and (H) we find: 

Ff = F°-(F°-F,°)(l-Pie-), (8) 

which shows that the overall conversion effect, both in the star and in the Earth, leads to 
suppression of the flux at low energies {E < Ec) and to its enhancement at high energies 
{E > Ec) ■ The effect is proportional to the difference of the original fluxes. 

Let us consider the properties of the Earth regeneration effect in some detail. The 
formula (H) can be written as: 

Ff = F° - (F° - F,° ) (sin^ 6 - A, sin^ 0) , (9) 

where 

Ap = sin 26^ sin {29 - 29^) , = — (10) 

'"m 

are the depth and the phase of oscillations. In contrast with the conversion in the star, 
the Earth matter effect has strong energy dependence. The depth of oscillations, Ap, 
equals zero at very low energies, where 9m — 9, and at high energies, where 9m — 7r/2. 
The depth reaches the maximum 

^;r"'' = sin^ ^ (11) 

at 9m = 9/2, which corresponds to the equality 

r = 1 • (12) 

Notice that the depth of oscillations increases with 9. At the condition of maximal depth, 
eq. (0), also the conversion probability Pig reaches the maximum: 

Ple-=1, (13) 

provided that = 7r/2 + kn with k = 0, 1, 2, . . . . For E -^ ot E -^ oo, we have 
Pie — >■ cos^ 9. Thus, Ap as well as Pig have a resonance character with eq. (12) being the 



resonance condition. Notice that the condition (|r^) does not depend on the mixing angle, 
in contrast with the resonance condition for the pure flavour case. 
From eq. (|12|) we flnd the resonance energy 

^ ~ 2V2G^e ' 



The interval of energies E^ -^ E^ around ( [14|) for which Ap > y4™"^'/2 = sin^^/2 is 
determined by: 

-i = cos 2^ + 2 ± a/(cos2^ + 1)(cos2^ + 3) . (15) 

Thus Er and E± give the interval of strong oscillation effect in the Earth. 

The peak in the Pig probability exists in the antineutrino channel even if the usual 
flavour resonance is in the neutrino channel. With increase of 6 and its shift to the dark 
side, 9 ^ 7r/4, the peak gets narrower and at ^ ^ 7r/2 it converges to the usual flavour 
resonance peak. 

The oscillation length l^ decreases with E/Am'^; it also decreases with cos 2^. 

At high energies, due to the exponential decrease of the fluxes with E and the difference 
of temperatures of the z/g and z/^ spectra, we get F?, ^ Eg, and therefore, according to 
eq. (I), the flux at the detector equals: 

Fi^^F°(l-Pie-). (16) 

Then, at the resonance, E ~ Er, eq. (^^, the Ue flux can be suppressed completely: 
FP ~ 0, if the phase of oscillations equals a semi-integer fraction of vr. This suppression 
does not depend on the value of the mixing angle 6. In what follows we will show that 
this suppression could indeed be realized in the case of neutrinos from SN1987A. 

3 Neutrino signals in K2 and 1MB 

Let us consider the difference of signals at K2 and 1MB which can be produced by oscil- 
lations in the Earth. This difference is related to the distances travelled by the neutrinos 
in the Earth: djMB — 8535 Km for 1MB, (1x2 — 4363 Km for K2, and to the average den- 
sities piMB — 4.5 g ■ cm~^, pk2 — 3.5 g ■ cm~'^ along the trajectories. As a consequence, 
both the depths and the phases of oscillations at K2 and 1MB are different. 

Let us show that the features (i)-(iii) indicated in the introduction can be explained by 
oscillations in the matter of the Earth. Moreover, the explanation implies certain values 
of Arri^ and sin^ 29. 

To reproduce the characteristics described in (i)-(iii) we require that: 

(1) The phase of oscillations at 1MB detector at ii^ ~ 38 ^ 42 MeV equals 

0/mb(4O) ^ ^^^ = fcvr , A; = 1,2,3,... . (17) 

I'm 

Under this condition the oscillations in the matter of the Earth do not suppress the signal. 

(2) The phase of oscillations at 1MB at ii^ ~ 50 ^ 60 MeV is a semi-integer of vr: 

0/MB(6O)=7r(^- + fc) , A; = 0,1,2,.. , (18) 
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so that in this range of energy one expects maximal suppression effect. It is easy to check 
that, in the range of parameters of interest, the conditions (0) and ([T8| ) are satisfied 
simultaneously with good precision. 

(3) The phase of oscillations at K2 at E ^ 38 ^ 42 MeV is 

0^2(40) = ^ = TT (^^ + fc) , A; = 0, 1, 2, ... , (19) 

so that the Earth matter effect produces maximal suppression of the K2 signal in the 
range E ~ 38 ^ 42 MeV. 

(4) The Earth matter effect is maximal, Ap ~ A™"^', at 1MB at the energies i? ^ 50 -^ 60 
MeV, that is: 

E^j^^ ~ 50 ^ 60 MeV . (20) 



In fig. we show the conditions (p!7|), (|T9|) and ( pO] ) in the IS.vp? — cos 20 plane. As 
follows from the figure, there are bands in which the requirements (|l^) and ([T9|) are 
satisfied simultaneously. They correspond to 4>imb — '^4'K2 = 37r, 57r, Tvr, ... . The phase 
increases with Am^. Notice that the requirements (p!^)-([T9|) are satisfied in the whole 
relevant range of cos 29 if 0/mb equals odd multiples of vr. 

The condition of maximal effect, eq. (pOD, gives Arn^ = (1.7 -r- 2.1) ■ 10^^ eV^. Large 
Earth matter effect, e.g. Ap ^ O.TA^"-^, is realized in much wider interval, whose borders 
depend on cos26' (the upper border is represented by the dashed line in fig. |I]). 

As follows from fig. |I]the explanation of the properties (i)-(iii) implies the oscillation 
parameters to be in the regions (piMB = 37r for all the relevant values of 9, and (J)imb = 57r 
for cos 26' ^ 0.2. The central values of these bands are described, approximatively, by the 
following lines: 

Am^ ^ 3.3 ■ 10^^ eV^ [1 - 0.35 cos 20] , 

Am^ ^5.6- 10"^ eV^ [1-0. 18 cos 20] . (21) 



For cos 20 ^ 0.1 the values (pTj) are well inside the 99% CL. allowed region of the LMA 
solution (dotted-dashed contour in fig. |1], from ref. [0). Moreover, the best-fit point, 
cos 20 ~ 0.5, Am^ ~ 5.2 • 10"^ eV^ |], lies in the band 4>imb — Stt (see eq. (pT))). 

The fig. |I| shows also the 99% C.L. exclusion curve obtained in ref. [^ (dotted line), 
corresponding to the upper bound on the permutation parameter p = 0.35. One can see 
that the region cos 20 ^ 0.25 is excluded; this region was obtained taking Tg ~ 4.5 MeV 
and T^i ~ 7.4 MeV; it becomes narrower for the lower values of the temperatures we 
consider in this paper. 

^This result is given by the global two- neutrinos fit of the solar neutrino data including the total rates 
and the day and night spectra at Super-Kamiokande ||lq] . 



The expected distributions of events are rather sensitive to variations of Am^. For 
instance, a change of Am^ by ~ 20% in the (piMB = Svr band will lead to the prediction 
of strong suppression of the number of events in the interval E' ^ 35 ^ 40 MeV of the 
1MB spectrum and enhancement of the signal at -E ~ 50 MeV, in contradiction with 
observations. 



> 
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Figure 1: Bands of equal phases 0/mb(4O) = krc (dotted regions) and 0i^2(4O) = 
TT (1/2 + k) (dashed regions) in the Ar7i^-cos2^ plane. The widths of the bands are de- 
termined by the requirement that the conditions (0) and (|17|) are satisfied in the energy 
interval ii^ = 38-i-42 MeV. The region below the dashed line represents the band of strong 
Earth matter effect, where Ap ^ O.TAZ^"'^ (see eq. (pO])). For comparison we show the 
99% C.L. allowed region of the LMA solution of the solar neutrino deficit (dotted-dashed 
contour, from W^) where the best fit point is marked by an asterisk. The dotted line 
represents the 99% C.L. exclusion curve from fig. 3a of ref. |jl2 . 



Let us consider the expected spectra of events at K2 and 1MB. The original instanta- 
neous z/e and z/^/ fluxes can be described by "pinched" Fermi-Dirac spectra with pinching 
parameter 77 ~ 2 [^]. However, once integrated over arrival time intervals At of several 
seconds, the spectra are well approximated by the ordinary Fermi-Dirac form, i.e. with 
r] = 0, as an effect of the decay with time of the neutrino luminosities and temperatures. 
This description is good if At is not larger than the typical decay time r of the temper- 
atures and luminosities. For At ^ r the integrated spectra can not be approximated by 



8 



a Fermi-Dirac form. This is the case of the SN1987A data, which show a rapid coohng 
of the neutrino spectra with decay time comparable or smaller than the duration of the 
burst, r ^ At^ ~ 13 s. Therefore, we have divided the whole time interval of observations 
in two bins, ti = 0^6.5 s and ^2 = 6.5-^13 s, and described the integrated fluxes over each 
bin by Fermi-Dirac spectra with different temperatures (T(ti) > T{t2)) and luminosities 
(L(ti) > L{t2)). We mark that the values of T and L we give should be considered as 
effective (time averaged) quantities. 

In figures ^^ we show the expected spectra of events at K2 and 1MB in the first time 
bin (t < 6.5 s) for two sets of parameters from the preferable regions ([21|). They have 
been obtained taking Tg = 3.5 MeV, T^/ = 7 MeV and 77^/ = T^g = 0. We assumed equal 
integrated luminosities in z/g and v^r. Lg = L^' = 3 ■ 10^^ ergs. The energy thresholds and 
detection efficiencies have been taken into account. 

According to the figures, the predicted spectra with oscillations (solid lines) fit better 
the observed distributions of events at K2 and 1MB. As expected, the spectra without os- 
cillations (dashed lines) can not describe the concentration of 1MB events at ii^ ~ 30 -h 40 
MeV and the absence of an excess of K2 events in this range. The conversion in the star 
only leads to an appearance of high energy tails (short dashed lines) which contradict 
the observations of both K2 and 1MB. The oscillations in the matter of the Earth (solid 
lines) suppress the tails above 45 MeV in 1MB and above 35 MeV in K2 |. Notice that 
the figures have illustrative character and do not correspond to the optimal (best fitted) 
choice of parameters of the original spectra. 

To quantify the improvement of the fit due to conversion inside the star and the Earth 
matter effect we have performed the likelihood analysis following, in general, the prescrip- 
tions of ref . |l^. There are, however, some differences: 



1). We have divided the whole energy range of the detected events into several bins. 
The size Ai? of each bin has been chosen according to the experimental errors e in the 
measurements of energy: /\E ~ 2e. We used three bins with Ai? = 10 MeV for K2 and 
two bins with AE = 15 MeV for 1MB (see histograms in the figs. ^|]). 

2). We computed the realization probability P of a given theoretical prediction using the 
Poissonian distribution: 

P = U ^e- , (22) 

where Ni is the number of observed events in the i-th bin and n, = ni{Te, Tp_i, Lg, Lp^/, Am^, 
6) is the corresponding number of expected events which depends on the parameters of 
the original neutrino spectra and oscillation parameters. We find the maximum of P and 



the contours of given confidence levels according to ref. ||T3 



^Notice that in our analysis the maximal suppression, due to oscillations in the Earth, of the antineu- 
trino signal in 1MB occurs at energies above 50 MeV. For such energies the efficiency of the 1MB detector 
was high {Eimb ^ 0.7). 



Notice that, in contrast with JlSl we use a discrete (binned) expression of the prob- 



abihty P and we calculate rii by integrating the corresponding predicted function n{E) 
over the i-th energy bin. By this procedure we take into account the experimental errors 



on energies with no need to introduce any energy resolution function as it is done in [ p!3[] . 

3). We performed separate analyses of the K2 and 1MB data collected in the two time 
bins ti and ^2; no events are discarded in this procedure^. 

In fig. H we show the results of the likelihood analysis of the 1MB and K2 data of the 
first time bin in the Tg — Lg plane in absence of mixing, ^ = (upper panel). As follows 
from the figure, K2 and 1MB imply substantially different temperatures of the original Ue 
spectrum: Tk2 ^ 2.8 MeV and Tjmb ^ 4.2 MeV. There is no overlap of the 68% C.L. 
regions. 

The lower panel of fig. ^ shows the analogous plot in presence of oscillations in the star 
and in the matter of the Earth. We have taken T^z/Tg = 1.8, L^//Lg = 1, rje = rj/ii = 
with oscillation parameters from the preferable bands (see fig. p. With oscillations the 
best fit parameters inferred from 1MB and K2 become closer: Tk2 — 1-96 MeV and 
TiMB ^ 2.75 MeV. We get AT/T = 2{Tjmb - Tk2)/{Timb + Tki) ^ 0.33 and 0.41 with 
and without oscillations respectively. Overlap of the 68% C.L. regions appears. The 
combined fit gives Tg ~ 2.77 MeV and Lg ~ 4.4 ■ 10^^ ergs, which are in good agreement 



with recent calculations ^7j. Thus, the likelihood analysis with oscillations shows some 
improvement of the fit and does not imply too low temperatures and high luminosities. 

In the second time bin K2 has detected 3 events with low energies, close to the thresh- 
old, whereas 1MB has no signals. This can be easily interpreted in the assumption of 
lower temperatures of the Ue and i/^' original spectra and smaller difference of T^/ and Tg. 
Taking for instance Tg ~ T^/ ~ 1.8 MeV and Lg ~ Lp_/ ~ 3 ■ 10^^ ergs we predict ~ 2 
events in K2 with energy ii^ ~ 5 -^ 15 MeV and less than 0.5 events at ii^ > 15 MeV. With 
the same values of the parameters we get less than 0.2 events in 1MB, that is, with high 
probability all the events are below the 1MB threshold. 

For neutrino parameters from the LOW solution: Am"^ = (0.3^2) ■ 10^^ eV^, sin^ 26 ^ 
0.95, we get Iq/Iu ^ 10~^, so that the mixing angle in matter is suppressed: sm29m — 
sin2^(/o//i/). For the depth of oscillations in the Earth matter we find from eq. (pi]|): 

(23) 

that is, Ap <^ 1. Thus, the effect of oscillations in the Earth is negligibly small. Now 
Pig ~ cos^ 6 and the i?e flux at Earth is determined by the conversion inside the star only 
(see eq. (^ and the short dashed lines in the figures 0-^. No improvement of the fit can 
be obtained in this case. 




^The event number 6 in K2 has been considered as due to background, following ref. ||l^; however 
our conclusions remain unchanged when this event is considered. 



10 



4 Implications for the neutrino mass spectrum 

The fact of observation of the Earth matter effect in the z/g channel by itself has important 
implications for the neutrino mass spectrum and mixing (see also |l( 



1). As we have shown, a significant effect is possible for oscillation parameters in the 
LMA region only: thus, the observation of the effect selects the LMA solution of the solar 
neutrino problem. In the cases of SMA or LOW solutions the Earth matter effect would 
be practically unobservable. 

2). Important conclusions can be drawn on the type of mass hierarchy and the mixing 
element Ues in the three neutrino context. Let us consider the three neutrino scheme which 
explains also the atmospheric neutrino problem via z/^ ^-» Uj. oscillations. In this scheme 
the third mass eigenstate, z^s, is isolated from ui and z/2, which are responsible for the solar 
neutrino oscillations, by the mass gap Am^s — Am^^^ ~ 3 ■ 10~^ eV^. The state u^ has 
nearly maximal mixture of i^^ and Ut-. The admixture of the z/g in this state, described by 
the matrix element |f/e3p, is small being restricted by the reactor experiments CHOOZ 
ig and Palo Verde [|l|: |f/e3p ^ 0.02 ^ 0.05. 



Observation of the matter effect means that (i) either the mass hierarchy is normal: 
Am23 = rn^ — m\ > (the third state is the heaviest one) or, (ii) if the hierarchy is 
inverted, Aml^ < 0, the z/g — z^s mixing is very small: |t/e3p -C 10~^. 

Indeed, in the first case the high density resonance associated with Am23 (pr ~ Aml^) 
is in the neutrino channel. It does not influence the conversion of antineutrinos and the 
situation coincides with the one considered above. 

If, however, the hierarchy is inverted the high density resonance is in the antineu- 
trino channel. The pattern of conversion then strongly depends on the adiabaticity in 
this resonance. If |?7e3p ~ 10^^ the level crossing is adiabatic, so that in the star the 



conversions 



z>; , z7; -^ z7' /z7e , and u' -^ v'Jv, (24) 



occur. Thus, the z/g flux at Earth will be composed of converted z/^ and v'^ fluxes. Since 

the original v' and z/^ fluxes are identical, no matter effect should be seen in the z/g flux 

TE| , pm . Moreover, the z/g flux will have the hard spectrum of the original v' and v'^ fluxes. 



which is disfavoured by observations [^ . 

If |f/e3p ^ 10^^, the adiabaticity in the high density resonance is broken. The effect 
of this resonance can be neglected and the pattern of the three neutrinos conversion is 
reduced to the two neutrino case discussed in the paper. 

If |f/e3p — 10~^-i-10~^, one can observe an intermediate situation: harder z/g spectrum 
and partially suppressed Earth matter effect. 

If |f/e3p ~ 10^^ (which looks rather plausible in view of the large mixings between 
the first and the second as well as the second and the third generations) we come to 
the conclusion that the observation of the Earth matter effect in the z/g channel implies 
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a scheme with normal mass hierarchy and the LMA solution of the solar neutrino problem. 

3). Let us consider the Earth matter effect for supernova neutrinos in the case of i/g — Ug 
mixing. As for active neutrinos, for Uf. — i^s ^ significant effect is possible for large mixing 
angles (notice, however, that large mixing z/g —>■ Ug conversion gives poor fit of the solar 
neutrino data). 

We assume that no sterile antineutrino fluxes are generated in the central part of the 
star^ In this case we get the result for the Ue flux at the detector from eq. (||) putting 
F',' = 0: 

FF = Pi-eFl . (25) 

Here Pfg is the probability of vi — > Ve oscillations in the matter of the Earth for mixing 
with sterile neutrino. 

The difference of potentials for the z/g — Vs system, Ves = \/2Gf [fie — nn/2) ^ V^e^/2, 
is about two times smaller than the one of the Ve and z/^ species, K^. In the last equality 
we considered that the medium is almost isotopically neutral in the mantle of the Earth. 
Correspondingly, the refraction length is two times larger than in the active-active case, 
so that, according to eq. ([T2|) the range of large matter effects is shifted to smaller Am? 
by a factor 2. 

From eq. ([25| ) we see that, in contrast with z/g — v' mixing, the observed events in the 
detector are due to the original Ve flux only, and no hard tail appears in the spectrum. As 
a consequence, no improvement in the fit of K2 and 1MB data is obtained with respect to 
the ^ = case. Moreover, it is possible to check that a larger Vf. luminosity is required. 

The discussion of the sections Bfl can be immediately generalized to the case in which 



the electron neutrino is mixed with both an active neutrino v^ and a sterile one, Vg. The 
problem can be reduced to Ue — ^x mixing, with v^ = Vl — ^^'/i + V^^'s! the pure active 
and pure sterile mixings correspond to 5 = and 5 = 1 respectively. Now the z/g flux in 
the detector is given by eq. (||) with F° = (1 — S)F^,. Thus, with respect to the pure 
active case, one gets a reduction by a factor ~ (1 — 5) of the hard part of the detected 
spectrum. For 6 ^ 0.5 still one gets a good fit of the K2 and 1MB signals. 

5 Conclusions 

We have studied properties of the oscillation effects in the matter of the Earth on neu- 
trino signals from supernovae. We show that certain features of the neutrino signal from 
SN1987A detected by K2 and 1MB detectors (difference of spectra, absence of events 
above i^ ~ 40 MeV) can be explained by oscillations in the Earth and different positions 
of the detectors. The consistency of the z/g spectra implied by K2 and 1MB data improves 
when oscillations are taken into account. 



■*A sterile state Vs can be produced, however, by conversion of active neutrinos at high densities in 
the star. This happens in some schemes with four neutrinos which explain the LSND result. 
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The existence of the Earth matter effects on the neutrino burst from supernovae has 
crucial imphcations for the neutrinos mass spectrum: 

(1) The oscillation parameters Am'^ and sin^ 26 should be in the range of LMA solution of 
the solar neutrino problem. Moreover, the values of Am^ and sin^ 29 should lie in rather 
narrow bands. 

(2) The mass hierarchy (in the three neutrino context) should be normal, if |t/e3p ~ 10^'^ 
and it can be inverted provided that |t/e3p -C 10^^. 

It should be stressed that, in view of the low statistics of SN1987A neutrino signals, 
our conclusions have an indicative character, as far as the interpretation of the K2 and 
1MB data is concerned. If the LMA solution will be identified in future experiments, 
the inclusion of the effects of oscillations in the matter of the Earth in the analysis of 
SN1987A data will be unavoidable, and the interpretation of SN1987A signals given in 
this paper will be confirmed. The analysis performed in this paper shows that searches 
for oscillation effects in the matter of the Earth in future observations of galactic super- 
novae will give important information on the neutrino mass spectrum, independently of 
supernova models. 
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Note added 



In the recent paper [^ a new analysis of the SN1987A signal is performed, based on the 
Kolmogorov test of the K2 data only, with the conclusion that the LMA solution of the 
solar neutrino problem is excluded at high confidence level. However the Earth matter 
effect, which we have shown here to be crucial, is not taken into account. 
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Figure 2: The predicted spectra of v^ events at Kamiokande-2 (upper panel) and 1MB 
(lower panel). We show: the original spectra without oscillation effects (long dashed 
lines), spectra with conversion in the star only (dashed lines), spectra with oscillation 
effects both in the star and in the Earth (solid lines). We used the following set of 
oscillation parameters and characteristics of the original z/g and v' spectra: Tg = 3.5 

0.5, 



MeV, Lg = 3 ■ 10^^ ergs, r], 



0, T^' 



7 MeV, L^/ = 3 ■ 10^^ ergs, v = 0> cos 2^ 



Am^ = 2.75 ■ 10 ^ eV^. The histograms show the observed distributions of events during 
the first 6.5 seconds. 
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Figure 3: The same as fig. with cos 20 = 0.2, Am^ = 3 ■ 10"^ eV^ 
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Figure 4: Best fit points and contours of equal 68, 90, 95.4% likelihood for K2 and 1MB 

in the Tg — Lg plane. The upper panel shows the result of the separate fits of K2 and 

1MB data without oscillation effects. The lower panel represents a similar fit in presence 

of oscillations. The following values for the oscillation parameters and characteristics of 

original spectra have been used: T^z/Tg = 1.8, Lp//Lg = 1, rjg = Vp.' = ^^^ cos 2^ = 0.2, 

Am2 = 3 ■ 10-5 eV2. 
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